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Abstract. Motivated by recent experiments on colloidal systems with competing attractive and repulsive 
interactions, we simulate a two-dimensional system of colloids with competing interactions that can un- 
dergo fragmentation. In the absence of any other confining potential, the colloids can form stable clusters 
depending on the strength of the short range attractive term. By suddenly changing the strength of one 
of the interaction terms we find a rich variety of fragmentation behavior which is affected by the existence 
of "magic" cluster numbers. Such soft matter systems can be used to construct artificial nuclei. 

PACS. 82.70.Dd Colloids - 45.70.Qj Pattern formation - 05.70.Ln Nonequilibrium and irreversible ther- 
modynamics 



The issue of pattern stability has attracted consider- 
able interest in a wide range of different areas. A com- 
monly encountered mechanism for pattern formation is 
the competition between short-range attraction and long- 
range repulsion, which occurs in many systems including 
magnetic films, Langmuir monolayers, diblock copolymers 
[1], water-oil mixtures [2], and two-dimensional electron 
systems [3]. Similar types of pattern formation are ex- 
pected to occur in dense nuclear matter, where they have 
been termed "pasta" phases. These phases may occur un- 
der conditions of stellar collapse [4] or in the crust of a 
neutron star [5] . Models of classically interacting particles 
with competing attractive and repulsive interactions have 
been used to study a variety of topics in nuclear mat- 
ter. For example, if two nuclei collide, they temporarily 
form a high density state which then expands rapidly and 
can fracture in a complex manner [BJ. Nuclear fragmen- 
tation processes in both two-dimensional [7] and three- 
dimensional [8] geometries cause the nucleus to be decom- 
posed into clusters of all sizes through mechanisms ranging 
from evaporation to multifragmentation. 

Models of small numbers of interacting particles, out- 
side the thermodynamic limit, can be useful in under- 
standing the structure and stability of naturally occurring 
clusters such as nuclei. Repulsively interacting particles 
can be forced to form clusters by confining the particles 
in a trap potential. A wide variety of systems exhibit this 
type of behavior, including electrons in quantum dots [9] 
or on the surface of liquid helium [10] , vortices in superflu- 
ids colloidal particles in circular traps [12], confined 
ferromagnetic particles [13] , and charged dust particles in 
plasma traps [H] . The resulting structures have been stud- 
ied extensively for a range of interaction types and trap 



types [15] . and the details are sensitive to the nature of 
the trapping potential. 

There has recently been growing activity in studying 
colloidal systems with competing long range repulsive and 
short attractive interactions, where various types of clus- 
ters are observed to form |16|17|18|19|20| . For instance, 
clusters and "artificial atoms" were created in an experi- 
mental system of rotating disks with no confining poten- 
tial [21] . Other current experiments use combinations of 
entropic attractive forces with additional repulsive inter- 
actions to form clusters and study their decay [22], and 
long-range repulsion between colloids has been achieved in 
nonpolar solvents [23] and water-oil emulsions [24] . Analo- 
gies have been made between colloidal clusters and nuclear 
matter, leading to proposals for mimicking nuclear matter 
phenomena with a colloidal system [20] . Although there 
have been several colloidal studies of equilibrium cluster 
shapes that have identified magic numbers, to our knowl- 
edge the fragmentation of colloidal clusters has not been 
considered previously. 

Here, we form clusters out of colloids that have a com- 
peting short-range attraction and a long-range Coulomb 
repulsion. We identify the stability line as a function of 
the strength of the short-range attraction for clusters of 
different sizes. When the attraction falls above this line, 
clusters exist in the absence of a confining potential. We 
then study different mechanisms for destabilizing the clus- 
ters, including thermal excitation and a sudden quench of 
the attractive term to a lower value. The latter process is 
inspired by the sudden expansion of an energized nucleus. 
The microscopic details of the decays can be accessed di- 
rectly through our simulations. We find numerous decay 
modes depending on the cluster size and the quench depth. 
This system can serve as a model for understanding the 
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Fig. 2. (a) Dots: Critical value B c for stable cluster formation 
versus cluster size N. Dotted line: Fit to B c = 4.2 + 0.5 ln(iV). 
(b) dB c /dN for the same data. Negative values of dB c /dN 
indicate clusters that are more stable than neighboring cluster 
sizes. 
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Fig. 1. Structure of stable clusters for N=(a) 5, (b) 6, (c) 7, 
(d) 8, (e) 9, (f) 10, (g) 11, (h) 12, (i) 13, (j) 14, (k) 15, (1) 16, 
(m) 17, (n) 18, (o) 19, (p) 20, (q) 22, (r) 24, (s) 26, (t) 30. 

complex interactions between charged fragments in a clus- 
ter that is breaking apart. 

We simulate N interacting colloids in a system with 
open boundary conditions. No confining force is applied 
to the colloids, which either remain in a cluster due to the 
colloid-colloid interactions when placed in a stable config- 
uration, or scatter under Coulomb repulsion when placed 
in an unstable configuration. The overdamped equation of 
motion for colloid i is r\dvij dt = f s ; = % + f T > where 

rj is a phenomenological damping term. The colloid-colloid 
interaction fy = — VJ7(r.y)fy for colloids separated by 
Vij = Ti — Yj consists of a long-range Coulomb repulsion 
and a short-range exponential attraction: 

U(r) = l/r - Bexp(-Kr). (1) 

At small and large r the repulsive Coulomb term domi- 
nates. The attractive interaction can be varied using the 
inverse screening length n and the parameter B. Here we 
fix k = 1. As we will demonstrate, clusters are stabilized 
above a critical value B c which is a function of N . The 
temperature is modeled as Langevin random kicks with 
the properties (f T (t)) = and {f \t)f \t')) = 2 V k B TS{t~ 
t'). To initialize the system, we set B > B c and allow the 
colloid positions to relax into a stationary state in the 
absence of temperature. 

Due to the attractive component of the colloid-colloid 
interaction, stable clusters of colloids form even in the 
absence of a confining potential. We illustrate the stable 
clusters obtained for representative values of N ranging 
from N = 5 to N = 30 in Fig. [TJ In general, the colloids 
tend to form a small portion of a triangular lattice when- 
ever possible, not unlike the hypothetical triangular struc- 
ture assumed in the nuclear lattice model [25]. Since the 
underlying substrate is flat, geometrically necessary dislo- 
cations do not arise. Certain configurations which are the 



most lattice-like are also the most stable. These "magic" 
clusters include N — 7, 12, and 19. Note that the magic 
cluster N — 12 forms a nonhexagonal segment of a tri- 
angular lattice, so that it is highly stable even though 
it does not meet the proposed magic cluster criterion of 
N — 1 + 3p(p — 1) for integer p [55]. The strong tendency 
toward triangular ordering is in contrast to systems of con- 
fined repulsively interacting particles, which tend to form 
ringlike structures instead. The smaller clusters shown in 
Fig. Q] resemble structures obtained in the cluster phase 
of either a periodic system |27] or one confined to a trap 
|28j . In these systems, interactions with neighboring clus- 
ters alter the shapes of larger clusters and can distort the 
triangular ordering. Thus, the configurations we observe 
here for larger unconfined clumps differ from those found 
in confined or periodic systems. 

To measure the stability of the configurations at dif- 
ferent values of N, we perform a series of simulations in 
which B is initially set to a value which gives a stable con- 
figuration, and is then abruptly lowered to a new value. 
We can identify when B < B c by observing the cluster 
break into two or more pieces. In Fig. [2{a) we plot B c 
versus N through N — 80. Configurations with B > B c 
are stable, while those with B < B c break apart. We find 
that B c increases approximately logarithmically with N, 
but that there are noticeable fluctuations which are par- 
ticularly pronounced for small N. This is more easily seen 
in a plot of dB c /dN, shown in Fig. HJb). Magic clusters 
appear as pronounced dips in dB c /dN. This general be- 
havior is similar to that observed in artificial atoms, where 
a confining potential is present [53] ■ 

The structure of our clusters is predominantly trian- 
gular, but we can define the equivalent of outer and inner 
shells of colloids in the clusters. Unlike systems of particles 
confined in traps, our shells are not circular. The evolu- 
tion of the shell structure is shown in Fig.[3fa), where the 
number of colloids N s h e ii in each of the first four shells is 
plotted as a function of cluster size N. The second shell 
first appears at N — 7, while the third shell arises for 
N > 18 and the fourth shell for N > 29. There is some 
tendency for the shells to favor an occupancy of N s h e ii — 6, 
as seen by the presence of a step in N s h e ii at N s hell = 6 in 
Fig. [3ja) for all three of the shells. These results resemble 
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Fig. 3. (a) The number of colloids N a heii in the outermost 
shell (•), second shell (□), third shell (♦), and fourth shell (A) 
as a function of N. (b) Disintegration temperature T m as a 
function of B for B > B c and five different values of N. •: 
N = 3. □: N = 6. ♦: N = 9. A: N = 12. X: N = 15. 

measurements in dusty plasma systems [14] . We expect 
the shell structure shown here to be generally applicable 
to the type of colloid-colloid interaction we have assumed, 
and not to be parameter dependent. 

We next consider how to destabilize our cluster struc- 
tures so that they fragment. First we study thermal desta- 
bilization. Here, we fix B at a stable value B > B c , and 
simulate the system for a fixed number of time steps at 
each value of T, increasing T until the thermal fluctua- 
tions break the cluster apart. We perform a series of runs 
at each value of B to obtain an average melting temper- 
ature T m as a function of B. The results for five different 
cluster sizes are shown in Fig. Etb). In each case, T m in- 
creases with B since the attractive energy of the cluster is 
increasing. The rate of increase steepens as N increases. 
The destabilization is triggered by the escape of at least 
one of the particles from the cluster, and can thus be re- 
garded as a first passage process that becomes more likely 
as T increases. In our overdamped system, it is difficult for 
a single particle to accumulate enough energy to escape 
the cluster at low temperatures, unlike the evaporation 
process that can occur in a massive system. This tends to 
stabilize the clusters at low temperatures. 

The structures can also be destabilized by starting 
from a configuration with B > B c , and then suddenly low- 
ering B to a value B* < B c . This type of process could oc- 
cur after two nuclei collide and form a highly compressed 
but unstable structure, which subsequently fragments. In 
a colloidal system, sudden adjustments to the form of the 
interaction potential might be achieved through magnetic 
interactions, as was considered for purely repulsive par- 
ticles in Ref. [30 , or through electric interactions, as in 
Ref. [3T]. We begin with configurations prepared slightly 
above B c and perform a series of quenches at small but 
finite T to different values of B*. We measure 20 quenches 
for each value of B* in order to accumulate statistics on 
the resulting fragment configurations. The clusters break 
into smaller clusters which are stable at B*. As B c — B* 
increases, the number of fragment clusters increases and 
the average size of a fragment cluster drops. 

Different clusters show different types of decay modes. 
We notate these modes using (j 1 31 , j 2 02 , . . .jm m ) , where ji 
is the number of colloids in clusters of type i, rij i is the 
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Fig. 4. (a) Decay modes for N = 5. P is the probability 
of observing the mode, and B* is the quench depth. •: (I s ); 
□ : (l 3 ^ 1 ); A: (l\2 2 ); O: (2 1 ,3 1 ); O (l 2 ^ 1 ); x: (l 1 ^ 1 ). (b) 
Decay modes for N = 12. •: (l 12 ); □: (l 11 ^ 1 ); X: (1 10 ,2 2 ); ♦: 
(1 6 ,2 3 ); A: (2 2 ,3 1 ,5 1 ); +: (1 4 ,2\3 2 ); ■: (l 5 ^ 2 ^ 1 ); A: (1 3 ,3 3 ); 
O: (l 2 ^ 2 ^ 1 ); o: (l\3\4 2 ). (c) Major decay modes for N = 
15. •: (1 7 ,2 4 ); □: (l 6 ^ 3 ^ 1 ); A : (1 5 ,2 2 ,3 2 ); o: (1 4 ,2\3 3 ); M: 
(l 3 ^ 2 ^ 1 ); V : (l 2 ^ 1 ); +: {l' & ,3* ,5*); ♦: (l 1 ^ 1 ); X: 
(2\3\5 2 ); O (l 1 ^ 1 ^ 1 ); ■: (l 2 ^ 1 ^ 1 ); A: (3 1 ,5' J 1 ); ▼: 
(l 1 ^ 1 ,? 1 ); *: (l 2 ^ 1 ); <: (l\7 2 ); (2 1 ,6 1 ,7 1 ); >: (^S 1 ). 
(d) Average multiplicity (M) of the decay as a function of B* 
for N = 5 (•), 12 (□), and 15 (x). 



number of clusters of type i, and there are m different clus- 
ter types. We also measure the average multiplicity of the 
decay at a given value of B, (M(B)} — (YliLi n j m (&))■ 
We illustrate the probability P of each decay mode for 
N = 5 as a function of B* in Fig. [4ja). Here we observe 
all six possible decay modes, and find two dominant de- 
cay modes: a breakup into one 1-particle cluster and two 
2-particle clusters, denoted (1 ,2 2 ), or complete disinte- 
gration into five 1-particle clusters, denoted (l 5 ). For the 
magic cluster N — 7 (not shown) , we find only a single de- 
cay mode of complete disintegration, (l 7 ), unless B* is ex- 
tremely close to B c . The N = 7 state has such strong sym- 
metry that the particles are not able to form subclumps 
as they break away from the main clump. Fig. H£b) shows 
the decay modes for the magic cluster N = 12, illustrated 
in Fig. dja). The two dominant decay modes are both 
symmetric: a breakup into three 1-particle clusters and 
three 3-particle clusters (1 3 ,3 3 ), illustrated in Fig. [5jb), 
and complete disintegration, (l 12 ), shown in Fig.HJc). At 
B* values falling in the transition region between these 
two regimes for the case of a cluster held at finite temper- 
ature, numerous decay modes occur with no single mode 
dominating. In the case of an asymmetric, nonmagic clus- 
ter such as N — 15, plotted in Fig. BJc), a series of de- 
cay modes appears. Images of representative decays are 
shown in Fig.[5^d-f). For B* close to B c , the cluster breaks 
into two large asymmetric pieces, such as (G 1 ,!} 1 ) shown 
in Fig. [5jd) . As B* is lowered, smaller pieces begin to 
appear, as in the case of (2 1 ,6 1 ,7 1 ) in Fig. E^e). For B* 
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Fig. 5. Images of decay, (a) N = 12, initial configuration, (b) 
N = 12, B* = 5, (1 3 ,3 3 ). (c) N = 12, B* = 4.75, (l 12 ). (d) 
N = 15, B* = 5.5, (6\9 X ). (e) N = 15, B* = 5.35, (2 1 ,6 1 ,7 1 ). 
(f) N = 15, B* = 5.1, (1 4 ,2 4 ,3 X ). (g) N = 19, B* = 5.55, 
(l 3 ^ 3 ^ 2 ^ 1 ). (h) N = 24, B* = 5.7, (l^lO 1 ,^ 1 ). (i) N = 26, 
B* = 5.6, (1 3 ,2\3\9 2 ). 



considerably below B c , the cluster disintegrates into many 
small pieces, as illustrated in Fig. [3(f) for (l 4 ^ 4 ^ 1 ). Sim- 
ilar patterns occur for larger clusters. In Fig. |5Jg), the 
magic cluster N = 19 undergoes a nearly symmetric de- 
cay around a single central particle into the configuration 
(l 3 ^ 3 ^ 2 ^ 1 ). A decay into two large pieces combined with 
two small fragments is shown in Fig. [5jh) for A = 24. 
The relatively symmetric but elongated N = 26 cluster is 
illustrated breaking into two symmetric large fragments 
with N = 9 along with a handful of smaller fragments in 
Fig. [Hi)- 

The average multiplicity of the decay (M), or the av- 
erage number of fragments into which the cluster breaks, 
increases as B* decreases. This is shown in Fig. |U[d) for 
A~ = 5, 12, and 15. The presence of two dominant decay 
modes for A" = 5 and 12 appears as steps in (M) as a func- 
tion of B*. The asymmetric cluster N — 15 has a much 
more complex structure of decay modes, with none of the 
modes strongly dominant. This produces a smoother in- 
crease of (M) with decreasing B* . The clusters reach the 
point of complete disintegration into individual fragments 
((M) = A) at B* = 4.78 for N = 5, B* = 4.76 for 
N = 12, and B* = 4.62 for N = 15. 

In conclusion, motivated by recent experiments in col- 
loidal systems with competing interactions, we study a 
system that can form stable clusters in the absence of any 
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confining potential. We have shown that some cluster sizes 
are magic, and correspond to a defect-free segment of a tri- 
angular lattice. The strength of the short-range attractive 
term required to stabilize the clusters increases logarith- 
mically with the cluster size. Clusters can be melted by 
increasing the temperature, or destabilized by suddenly 
quenching the attractive portion of the interaction to a 
lower value. We find a variety of decay modes in the lat- 
ter case, with a structure that depends on the symmetry 
of the original cluster. These results suggest that colloidal 
systems of self-confined clusters where the interactions can 
be controlled may offer insights into the effects of interac- 
tions between charged fragments during decay processes 
such as multifragmentation for a larger class of systems 
including molecular clusters and nuclei. 

We thank W.K. Kegel for useful comments and bringing to 
our attention some of the recent experimental work on colloidal 
systems with competing interactions. This work was supported 
by the U.S. Department of Energy under Contract No. W-7405- 
ENG-36. 
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